MODELING A PACKED BED REACTOR UTILIZING THE SABATIER PROCESS 


Malay G. Shah, Anne J. Meier, Paul E. Hintze 
NASA, Kennedy Space Center, FL 32899 


ABSTRACT 


A numerical model is being developed using Python which characterizes the conversion and 
temperature profiles of a packed bed reactor (PBR) that utilizes the Sabatier process; the 
reaction produces methane and water from carbon dioxide and hydrogen. While the specific 
kinetics of the Sabatier reaction on the Ru/Al203 catalyst pellets are unknown, an empirical 
reaction rate equation? is used for the overall reaction. As this reaction is highly exothermic, 
proper thermal control is of the utmost importance to ensure maximum conversion and to 
avoid reactor runaway. It is therefore necessary to determine what wall temperature profile 
will ensure safe and efficient operation of the reactor. This wall temperature will be maintained 
by active thermal controls on the outer surface of the reactor. 


Two cylindrical PBRs are currently being tested experimentally and will be used for validation of 
the Python model. They are similar in design except one of them is larger and incorporates a 
preheat loop by feeding the reactant gas through a pipe along the center of the catalyst bed. 
The further complexity of adding a preheat pipe to the model to mimic the larger reactor is yet 
to be implemented and validated; preliminary validation is done using the smaller PBR with no 
reactant preheating. When mapping experimental values of the wall temperature from the 
smaller PBR into the Python model, a good approximation of the total conversion and 
temperature profile has been achieved. 


A separate CFD model incorporates more complex three-dimensional effects by including the 
solid catalyst pellets within the domain. The goal is to improve the Python model to the point 
where the results of other reactor geometry can be reasonably predicted relatively quickly 
when compared to the much more computationally expensive CFD approach. Once a reactor 
size is narrowed down using the Python approach, CFD will be used to generate a more 
thorough prediction of the reactor’s performance. 


NOMENCLATURE, ACRONYMS, ABBREVIATIONS 


CFD Computational fluid dynamics 
PBR Packed bed reactor 

Ru/Al203 Ruthenium on aluminum oxide 
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Project Background 


¢ What is the Sabatier process? 
CO, + 4H, - CH, + 2H2O 


kJ 
AH =—165.0™/, 
— Reactor utilizes Ru/Al,O, catalyst pellets 
¢ Very high selectivity for CH, production 
¢ Sample Reactor Design [1] 
— Preheat feed gas using excess heat 
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Preheat Loop 


Application: 

— Fuel on Mars! 
¢ Why is modeling important? 

— Estimate reactor hot spots 
Implement proper thermal control 
— Prevent catalyst degradation 
Predict scaled-up performance 


Modeling methodology: 
— Simplified 1D model written in Python 
¢ Very quick and gives good estimates 
¢ Good initial guess to further investigate using CFD 
— Traditional CFD 
¢ Longer to simulate, but may provide better reactor information 
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Model Information 


Reactor Dimensions 
— 0.43” diameter 
— 1.75” & 3.5” length 


— 0.125" x 0.125” Cylindrical pellets 
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¢ Assumptions 
— Steady-state 
— 1D porous media approximation 
— No intermediate reactions 
Catalyst effectiveness factor = 1 


« Reaction Rate 


— Ruthenium Catalyst: Empirical correlation [2] 


r = kg exp(—E,/RT) [pco,|"[pu, |” - ona] WPusol “Pol : /(T) 


— Nickel Catalyst: Hougen-Watson type [3] 
legkco. Ke | 
_ 14Xco,%y,Pco,PuH, (1-9) 
(1+ Kco,Pco, + Ky,Pu, ) 


1 Pcu,Pi,0 
Ke(T) Deo, Dit, 


B= 
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Model Information 


¢ Pressure Drop 
— Bed-to-particle diameter ratio is 3.4, so wall effects are not 
entirely negligible 
— Use Reichalt correlation instead of Ergun [4] 
¢ Heat Transfer 


— Thermal resistance method 
¢ Convection along inner wall 
¢ Effective radial conduction through bed 


— Many correlations! 
¢ Martin & Nilles: Convective heat transfer coefficient 
¢ Winterberg: Effective bed thermal conductivity 
¢ Zehner & Schlunder: Effective stagnant bed thermal conductivity 
1 1 1/Bi+3\ D, 
tae ne (= + 7) ae 
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Model Information 


¢ Differential Equations 
—- — Yk 
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Model Information 


Partial Pressure, 
pi 


Reaction Rate, r 


Internal 
Feed Line 


: Heat Transferred 
-----p} Between Reactor 


Heat of the : and Feed, q 
Reaction PP SET ee 


Effective Heat 
Transfer Coefficient, 
Vert 
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Thermal Conductivity 
Viscosity 
Density 
Average Velocity 
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Experimental Data 


Catalyst CO,/H 


mL SCCM A C E ee bed ene bed 
el el 
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Model Results 


bed volume | flow ee Be 
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| 23 | 40/200 | 225 | 295 | 314 | 295 | 258] 376 =| 293 
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Model Results 


bed volume | flow ee Be 


SCCM Cae bed ae bed i 
middle outlet 


| 23 | 120/600 | 93 | 282 | 366] 333 | 219] 523 | 337 
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Model Results 


bed volume | flow ee Be 


SCCM ere bed el bed mel 
ee —- 


| 2.3 | 160/800 | 85 | 298 | 416 | 395 _| 267 | 
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Model Results 
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Model Results 


bed volume | flow ee Bo 
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Model Results 


bed volume | flow ee Bo 


SCCM tae bed ee bed a 
middle outlet 


| 6.3 | 425/2125 | 253 | 467 | 530] 512 [475] 579 | 500 
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¢ Additional Experiments: 
— Larger reactor currently being tested 
— Radial profile to be captured 


¢ Python model: 
— Increase complexity 
¢ Radial profiles important! 
¢ Catalyst effectiveness 
— Refine boundary conditions 
— Validate against additional tests 


¢ CFD Model: 


— Catalyst surface reaction 
¢ May require more information on reaction kinetics 


TFAWS 2017 — August 21-25, 2017 16 


Preliminary CFD Model 


¢ Software: STAR-CCM+ 

¢ Catalyst pellets injected using Discrete Element Method 
(DEM) particles 
— Fluid volume gathered from subtracting resulting solid volume 

¢ Species Source Terms 
— Reaction modeled using field functions 


: \ + a> an ee ee e507 a 
DD See Soe 


Mole Fraction of CH4 
0.10000 0.15000 0.20000 0.25000 


0.0000 0.050000 


Temperature (C) 
461.69 510.98 609.56 
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